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Abstract 

r- ; 

£SJ , Non-lcptonic D meson decays may provide a reliable testbed for the multiquark interpretation of 

light scalar mesons. In this letter we consider D s decay and show that a 4-quark /o(980) meson could 
induce a decay pattern, which is forbidden for a qq constituent structure. Experimental tests to probe 
■ such possibilities are within reach in the near future. 

<N 

m I Preprint No. Romal-1449/2007 

Keywords Hadron Spectroscopy, Exotic Mesons, Heavy Meson decays 
^ I PACS 12.39. Mk, 12.40.-y,13.25.Jx 

Ph! 1 Introduction 

Recent observations of the three body non-leptonic decays of charmed mesons and of the J/ip strong 
decays in similar channels are giving important information on low-energy meson dynamics. We refer in 
particular to the irir and Kir mass distributions in D decays by E791 [I], BaBar [2], and to the similar 
distributions in J/tp — > (f)irir and <f>KK by BES [3]. 

In this paper, we analyse the S-wave amplitudes of D s decays: 



D s ^tt + (K + K-) (1) 
D s —* 7r + (KsKs) (2) 
D s -> TT+vr+vr^ (3) 

which can be observed by BaBar and Belle, following the study of the three Kaon decays of D° already 
published in [2]: 

D° K S K + K~ (4) 

In high statistics experiments, the S-wave amplitudes in ([T]) to (j4]) can be obtained from the Dalitz plot 
distributions [2j and are expected to be dominated by scalar meson exchange, fo and ao in particular. 

We point out that the observation of the /o(980) structure in the S-wave amplitudes of ([1]) and (|2|) 
offers a unique possibility to elucidate the valence quark composition of the light scalar mesons. The ratio 
of decay rates (|2|) to (pQ) is predicted to be close to 1/2 if /q(980) is a / = 0, qq state, while a tetraquark 
composition, /o = [sg][fg] (q = u,d) [1], could give a different result due to possible interference between 
1 = and 1=1 amplitudes. We show that, by fixing one single constant, the ratio of the two isospin 
couplings, the neutral K$Ks rate could be made almost to vanish in the accessible region. 
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The tetraquark nonet of scalar mesons is supposed to be completed by the /o(600)(=<r) [HE] and by 
the k(800) [7] resonances and we comment as well on the decays: 

D s vr/o(600) -h- 3tt (5) 
D s -» «if -»■ vrifif (6) 

Concerning decays ([!]), we explain the remarkable equality of the KK mass spectra for the charged 
and neutral combinations, K$K + and K~K + , as due to the rapid opening of the kaonic decay channels 
of /o and ao above the KK threshold, a mechanism pointed out long ago by Flatte et al. [5]. 

The results of the present work are summarized as follows. 

• The observation of a ratio much smaller than 1/2 for the S-wave rates ([2]) to (pQ) would provide a 
clear demonstration of the non-conventional nature of /o and ao mesons. 

• Once the ratio of the rates in ([2]) and (pQ) is given, one is led to a univocal prediction of the S-wave 
decay rates of (|3|) and of : 

D s -> ir°K + K s , ir°TT + r) (7) 
which we present in this paper, for the case of maximal interference. 

• Suppression of the a channel © is implied by the quark diagrams of Fig. [2] while the k band should 
in principle appear in the Dalitz plot Q, once the K* amplitude is subtracted. 

• We find that the experimental distributions of D° decays, Eq. (j4]), are consistent with the presence 
of both /o and Oq in fl3J), which is required by the four quark model and gives a consistent picture 
of the KK, S-wave production in D decays, see Fig. [5J 

We describe resonant amplitudes with a modified Breit-Wigner, according the the prescription of [5]. 
We use masses and couplings from [3], for the /o, and from [10] and [2] for the ao- 

For tetraquark states, an isospin violating fo-ao mixing could be induced by the u-d quark mass 
difference [S]. Data available thus far, however, do not show substantial mixing, see however [9], and 
we restrict to exact isospin for simplicity. Generalization of our formalism to the mixed case is straight- 
forward. High statistics data on decays (|3j) and (J7|) would provide a very sensitive determination of the 
mixing angle. 



2 Isospin structure of the weak amplitudes 

We consider the (bare) weak hamiltonian corresponding to charm-changing, parity-conserving, Cabibbo 
allowed decays (|H2p : 

H W = GC °^ 6C [(*7"c)(*v0 + (7m - 7m7s)] (8) 

Quark diagrams leading to final states with two quark pairs are reported in Fig. [H and include c- 
quark decay, (a), and annihilation, (b) to (d). Diagram (a) leads clearly to a KK state with 1 = and 
so does diagram (b), if we assume that the additional quark pair is produced from an isospin invariant, 
charge-conjugation symmetric sea. 

Diagrams (c) and (d) differ by the exchange u — > d in the weak current and d — > u in the sea quarks. 
Only axial weak currents contribute, due to the pseudoscalar nature of D s and to parity conservation. 
Therefore, for a symmetric sea, the amplitude takes a plus sign under these exchanges, leading again to 
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Figure 1: Quark diagrams leading to ir + plus one quark-antiquark pair; (a): quark decay. (b,c,d): quark annihilation. Only 
Axial x Axial currents contribute to diagrams (b) to (d), since the initial meson is pseudoscalar and the overall transition 
is parity conserving. 



an / = state for KK: 



A(D S 
A(D S 



n+K + K-) = A a + A b + A c 
= A a + A b + A d 



= A a + A b + A AA = Ai= 
A a + A b + A AA = Ai =0 



(9) 



Thus, taking into account the effect due to the difference of K + K and K°K° thresholds, the 
conventional qq picture of the scalar mesons implies uniquely: 



Rate(K s K s ) 
Rate( K+K~) 



1 



S—wave 



Rate(A°A°) 
Rate(K+K~) 



S—wave 



- x 0.83 



(/o,ao 



(10) 



The situation is different for the tetraquark structure, where we may expect that quark diagrams 
with three quark pairs in the final state dominate. The corresponding diagrams are reported in Fig. [2 
where we restrict to quark decay amplitudes for simplicity. In this case both axial-axial and vector- vector 
amplitudes contribute. The latter take a minus sign under the exchange u — > d. Therefore, we may have 
both / = and / = 1 for the KK: 



A(D S 
A(D S 



TT+K+K-] 
ir+K°K°) 



= A a + A b 
A a + A c = 



= A a + (A AA + A vv ) 
A a + (A AA - A vv ) 



= A I=0 + Ai=i 
Ai=o — Ai=\ 



where 



Ai = {(KK)j; out\Hw\D s 



111 



(12) 



A near cancellation is obtained if the axial-axial and vector-vector amplitudes are close to each other 
and the amplitude to produce the ir + from the weak current is negligible, a not unlikely situation since 
the latter amplitude vanishes with the pion mass in the free quark limit. In the following, we analyze the 
extreme case where negative interference in the K$K$ channel is maximal and (see SectH]): 



Rate(A 5 A 5 ) 
Rate(A+A") 



(/o,«o = [sq][sq] 



(13) 
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Figure 2: Quark-decay diagrams leading to ir + plus two quark-antiquark pairs. Contributions from Vector x Vector and 
Axial x Axial currents only, since the overall transition is parity conserving. 

3 Pole model of the S-waves 

The K + K~ mass distribution in ([T]) is dominated by the 0(1020). Once the dominant contribution is 
subtracted, the S-wave amplitude immediately below the 4> region should display the behavior associated 
with the rising Breit-Wigner of the /o(980) and ao(980) toward the peaks, which are however below the 
K K threshold. The <f> peak is, of course, absent in the K$K$ mass distribution, which should display 
the /o/ a o structure without any subtraction. 

As seen before, the conventional qq picture of the scalar mesons implies the KK system to be in pure 
1 = state. This needs not to be true if there are two independent isospin amplitudes for processes (pQ 
[2]) , as it is the case if /o and ao are tetraquark states [4] . 

We dominate the S-wave amplitudes in (jlip with the resonant processes: 

^^vr + / (980);7r + a (980); 

f /a ^K + K-;K°K° (14) 
Assuming exact isospin in the mass eigenstate^] we write: 



\f) = \S ) 
la) = I Si) 



V2 

[su][us] — [sd][ds] 



V2 

In an abbreviated notation, we can write the amplitudes of (pQ) and ([2j) according to: 



(15) 



A{D S ^it + KK) = (KK;out\ \^2\S I )BW I (s){S I \j \A) (16) 

and we introduce: 

g I = (KK;out\S I ) (17) 

M = (Si\A) (18) 



'we denote by e.g. [su] the fully antisymmetric combination of strange and up quarks [4] 
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with I = 0,1. The pole model is described by the diagram in Fig. [3j 




Figure 3: Pole diagram for Df, S-wave, decay amplitudes. A similar diagram has to be added for ao exchange. In the 
notation of Eq. Jl6j of text, A = (f\A) and g = (KK;out\f). 

BW/(s) and BW a (s) describe the line-shape of the resonances, which we take of the relativistic 
Breit-Wigner form, e.g.: 

BW > ~ s - Mj +iMTf{s) (19> 

Mf and Tf are the mass and total width. According to Ref. [5], we introduce an s-dependence in the 
widths, to take properly into account the opening of the KK threshold inside the fo or ao widths. 
The widths that appear in (|19|) are parameterized as follows. 

• The I =1 charged scalar states, that is the a^^SO), have a distinctive t]tt decay: 

0^(980) r ? vr ± (20) 

which, of course, is shared by the 1=0 component, S\. Conventionally, the corresponding rate Ti 
is written according to: 

ri = r(4(980) - n*±) = <4^# (2i) 



where p vn is the decay momentum (there would be no real need to introduce an s-dependence here, 
because the threshold of the final state is quite below the resonance position). In Ref. [10], T\ is 
parameterized with a coupling (gi)cB- 

( gi ) 2 CB = M a g 2 vn (22) 

We obtain from Table [TJ 

ri(M a )~69MeV (23) 

The 1=0 state, So, has a distinctive tttt decay. We denote by To(s) the 2ir width and follow the 
convention of Ref. O [5] : 

3 2 2 p n (s) 
2 V s 

where the factor 3/2 arises from the charge multiplicity of the final ir + ir~ and vr 7r° states and 
is the decay momentum, (also in this case the s-dependence is hardly needed). Using Table [2) 

T (M f ) ~ 236 MeV (25) 
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So,i have K + K and K°K° decays characterized by the couplings (fTHj) , We write, e.g. for fo: 

r(/ -> + r(/ -» = 

(26) 



9l 



2pch(s) 2p neu (g) 



We denote by p c h,neu an d Pol the decay momenta of K + K , K°K° and l^+i'f in the scalar particle 
decays, : 

Pch.neu = ~\Js - 4rn|T~ poi = ^=^A(s,m^ + ,m^ ) (27) 

and A(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2zx — 2zy, the usual triangular function. 
• We take go from Ref. [3]: 

dl = (9 2 k)bes - 694 MeV; (28) 
Refs. [10] and [2] give different determinations of g\\ 

2 _ [(92)cb} 2 _ j 54 MeV : Crystal Barrel 
91 2M a ~ \ 108 MeV : BaBar 1 j 

We take the BaBar value for consistency. 
Summarizing, we give the expressions of the widths for fo, a[j and as follows. 

3 2 2 P^( S ) 2 2 
2 

r a« = r i + 9l (Pch + Pncn) 

r a ±=r 1 + 2 5l — _. (so) 



Table 1: Parameters of the Oq , from [TO] and [2]. The subscript CB refers to the Crystal Barrel definition of the coupling, 
see Eq. (|22p . (gi.2)cB refer to the tjit and KK channels, respectively. 



M a (MeV) [10J 


( gi ) CB (MeV) HD1 


(92/9'ihB m 


(92)cb (MeV) 


999 ±2 


324±15 


1.03±0.14 


473±29±40 



Table 2: Parameters of the fo, from [3]. 



Af) (MeV) 


9i (MeV) 




965±8 (stat)±6 (syst) 


165±10±15 


4.21±0.25±0.21 
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The amplitude in (|16|) . where we use (|T9|) for the propagator of the unstable particle, corresponds to 
the Feynman diagram in Fig. [3] (plus the analogous one for ao exchange) . We obtain the following form 
for the differential distribution in s = m 2 ^: 

K K 



dT(D s -» ir+KK) 
ds 

1 - ^ (Ml - SM \ A{D , - , W ?m _ (31) 



128^ 3 M2 2Md s 1 v s " \/i KK 



where is given by (|z) 

With the same normalization of the amplitude, the Dalitz plot density is given by: 

dT(D -> abc) _ 1 , 2 



dsids 2 256vr 3 M 3 - 



A(£>->a&c)| (oi( )^ fll | 2 (32) 



£> (bc)" 

4 D s -> (7rif A^^g and D s -> (A+t?),^^ 

Above threshold, the KK widths grow rather quickly. One expects [5] the line shape to be dominated 
by the KK widths and therefore a quite similar behavior for /o and ao- 
Therefore, it is sufficient to tune the ratio: 

Q = ^ (33) 

to suppress almost completely the K°K° amplitude in the accessible region, independently from the fact 
that Ai=o and Ai=i are dominated by different (close lying) resonances. 

We analyze the decays on the basis of Eq. (f3Tj) . We take the values of the parameters for / from [3] 
and for from [10\ [2] , as summarized in Tables [1] and [2] . 

We report in Fig. HJ the theoretical curves for the phase space corrected probability, \F K ^-\ 2 Eq. ([3Tj) . 
for D s — ► 7r + K + K~ , D s — > ir + K°K° and D s — > ir°K + K° as functions of -y/i = m K j^. Curves are given 
for the central values of the parameters and for: 

Q = ^ = 1-5 (34) 

which minimizes the ratio of the rates ([2]) to ([I]), each integrated between threshold and 1.15 GeV. The 
suppression of the neutral channel is evident. 

The value in (I34p follows in fact from a very simple argument. Close to threshold we neglect the 
s — M 2 term as well as the kaon rates in the denominator of (1191) . so that: 



a «=Gfer o; ^ K TO - 41 (35) 

Using the numerical values, we see that the factors in parentheses are in the ratio of about 1.3 for 
(/ = 1)/(J = 0), so that the value (|34p leads to nearly cancel the K°K° amplitude at threshold, see 
Eq. <HH> 

With Q fixed, we can predict the S-wave rates of D s — > 37r and D s — > ?727r, ((7J), normalized to 
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Figure 4: Left panel. Phase-space corrected, S-wave probabilities, \F\ 2 see Eq. (|31[) . for the decay modes Dt — * ttKK 
as function of the invariant mass of the systems in parentheses. The value Q = Aq/Ai = 1.5 is assumed. Right panel. 
Distribution of the rates of D s — > 7v + (iT + iT~)h,i (red) and D s — » ■K°{ji + rj) (blue). The 7r + (7r ?7) distribution equals the 
tt°(it + t]) one. Absolute normalizations correspond to Ai = 1 in Eq. (|16[1 , For two pions we give only the high mass mmgh 
distribution. 



The 3-7T and r]2-rr decays are described each by two diagrams, differing by the exchange of 7t + tt + and 
7T + 7T°, respectively. We obtain the mass distribution by one-dimensional integration of the Dalitz-plot 
density, Eq. ([32]) . We give in Fig. HI right panel, the distribution for D s — > ir + (it + ir~) and D s — ► 7r°(7T + r/) 
as functions of the invariant mass of the particles indicated in parenthesis. The D s — > ir + (-K°r]) distribution 
is equal to the 7r°(7T + ry) one. 

In the limit of exact isospin (i.e. neglecting /o-ao mixing) the 3ir channel selects the 7 = scalars, 
/o(980) and /o(600). However the presence on one ss pair in the final state, Fig. [21 implies that the 
/o(600) should be suppressed as assumed in Fig. |H 

In conclusion, we report in Table 3 the ratios of the branching ratios of the D s modes considered thus 
far to the it + K + K~ mode: 

R(X) = ^ Br ' D -7 X ' , (36) 



Table 3: The numerical values of the ratios R(X), Eq. (|36[1 . for different values of Q 



Q 


R(tt + K°K°) 


R{^K + K") 


R{lT + 1T~-K + ) 


i?(? ? 7T U 7T + ) 


l 


0.23 


0.33 


1.12 


8.9 


1.5 


0.027 


0.23 


1.74 


6.15 


2 


0.046 


0.17 


2.3 


4.51 



5 D° K°(K + K-) S - W ave and D° K-(K+K ) S - W ave 

We consider these decays as arising from the quark process similar to those in Fig. [21 with an ss pair 
taken from the sea: 

C + Aspect -> S + (ud) + (SS) + -U spect (37) 

The final state mesons are formed from quark states by different diagrams with respect to those considered 
for D s decays. This indicates that the amplitudes of D s and D° are not related by symmetry, even in 
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Figure 5: Theory (for Q — 1.5) and data (boxes and triangles [2]) as functions of m K + K ~, blue dashed, and m K +^o, 
continuous red (curves normalized to the data). 

the exact flavor-SU3 limit. The weak hamiltonian (J8]) behaves under SU3 according to: 

H w = (sc){ud) = 6 15 (38) 

and there are altogether seven independent couplings for the S-wave trasition D — > Pseudoscalar +Scalar, 
with both S and P octet (eleven couplings, if we allow for singlets). The independent ways to couple the 
three quark pairs in Fig. [2]to the final qq octet meson are less (seven in all, including singlets), nonetheless 
there is no simple relations between the D s and D° couplings even in the quark diagram approximation. 

The K~ (K + K°) channel. The K~ is formed from the spectator u plus a strange quark taken either 
from the weak decay or from the sea. The amplitude is not related by symmetry to the D s amplitudes 
previously found. The system recoiling against the K~ is dominated by the a^. The dependence from 
yfs = m K x is well reproduced [2] with the parameters given in Table [H 

The K°(K + K~) channel. The K° is formed from the current d plus a strange quark taken either from 
the weak decay or from the sea. In both cases, a (susu) system is left, which should decay in K + K~ 
but not K°K°. We encounter here the same situation that we analyzed in Sect. [H We are unable to 
relate by symmetry the 1 = and 1 = 1 amplitudes of this decay, A' 1; to those of D s decays, eqs. ([18]) . 
However, in the pole approximation, the conspiracy between the /o and ao terms must be the same as 
in ([Ml), to cancel the K°K° rate. 

The K°K°K° channel. A clear prediction of the present scheme, is that this channel should be 
suppressed, if so is the D s — > -k + (K°K°). 

To compare with BaBar data [2j, we consider the phase space corrected distributions, \F K ^-\ 2 defined 
as in eq (j3ip . for the reactions ([4]). From the previous discussion, we derive the forms: 

A [D° - K-(K + K )] = C 91 BW a (s) 

A [ D o ^ k°(K+K-)] = C" [ gi BW a {s) + Qg BW f (s)} (39) 
with C and C" two independent constants and Q given by (I34h . 
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We report in Fig. [5] theory and data as functions of m K + K -, blue dashed, and m K +^o, continuous 
red (curves normalized to the data). 

The agreement is remarkable. In particular, it is to be noted the near coincidence of the K + K° and 
K + K~ distributions, which are dominated by a$ and by a superposition of /o and a[j, respectively. This, 
of course, is due to the validity of the argument of Flatte et al. [3], similarly to the D s case. 
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